arXiv:l504.04218vl [astro-ph.GA] 16 Apr 2015 


Mon. Not. R. Astron. Soc. 000, [T]-?? () 


Printed 17 April 2015 (MN L5 TbX style file v2.2) 


The Stellar Spectral Features of Nearby Galaxies in the 
Near-Infrared: Tracers of Thermally-Pulsing Asymptotic Giant 
Branch Stars? 

Rogerio Riffel^*, Rachel E. Mason^, Lucimara R Martins Alberto Rodriguez-Ardila^, 
Luis C. Ho^’^, Rogemar A. Riffel Paulina Lira Omaira Gonzalez Martin 
Daniel Ruschel-Dutra Almudena Alonso-Herrero^*^, Helene Llohic^\ 

Richard M. McDermid^’^^d6^ Cristina Ramos Almeida^’^^, Karun Thanjavur^^ and 
Claudia Winge^'^ 

1) Departarnento de Astronomia, Universidade Federal do Rio Grande do Sul. Av. Bento Gongalves 9500, Porto Alegre, RS. Brazil. 

2 ) Gemini Observatory, Northern Operations Center, 670 N. A ’ohoku Place, Hilo, HI 96720, USA. 

3) Laboratorio Nacional de Astroflsica - Riia dos Estados Unidos 154, Bairro das Nagoes. CEP 37504-364, Itajubd, MG, Brazil. 

4) Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China. 

5) Department of Astronomy, School of Physics, Peking University, Beijing 100871, China 

6) NAT - Universidade Cruzeiro do Sul, Rua Galvdo Bueno, 868, Sdo Paulo, SP, Brazil. 

7) Universidade Federal de Santa Maria, Departarnento de Fi'sica/CCNE, 97105-900, Santa Maria, RS, Brazil. 

8) Departarnento de Astronomia, Universidad de Chile, Casilla 36-D, Santiago, Chile. 

9) Institute de Astroflsica de Canarias, Calle Via Ldctea, s/n, E-38205, La Laguna, Tenerife, Spain. 

10) Institute de Fisica de Cantabria, CSIC-UC, 39005 Santander, Spain. 

11) University of the Pacific, Department of Physics, 3601 Pacific Avenue, Stockton, CA 95211, USA. 

12) Australian Astronomical Observatory, PO Box 296, Epping, NSW 1710, Australia. 

13) University of Victoria, Victoria, Canada. 

14) Gemini Obsen’atory, Southern Operations Center, c/o AURA, Casilla 603, La Serena, Chile. 

16) Department of Physics and Astronomy, Macquarie University, NSW 2109, Australia. 

17) Departarnento de Astroflsica, Universidad de La Laguna, E-38205, La Laguna, Tenerife, Spain. 


ABSTRACT 

We analyze the stellar absorption features in high signal-to-noise ratio near-infrared 
(NIR) spectra of the nuclear region of 12 nearby galaxies, mostly spirals. The features de¬ 
tected in some or all of the galaxies in this sample are the TiO (0.843 pm and 0.886 pm), VO 
(1.048 pm), CN (1.1 pm and 1.4 pm), H 2 O (1.4 pm and 1.9 pm) and CO (1.6 pm and 2.3 
pm) bands. The C 2 (1.17 pm and 1.76 pm) bands are generally weak or absent, although C 2 
(1.76 pm) may be weakly present in the mean galaxy spectrum. A deep feature near 0.93 pm, 
likely caused by CN, TiO and/or ZrO, is also detected in all objects. Fitting a combination 
of stellar spectra to the mean spectrum shows that the absorption features are produced by 
evolved stars: cool giants and supergiant stars in the early- or thermally-pulsing asymptotic 
giant branch (E-AGB or TP-AGB) phases. The high luminosity of TP-AGB stars, and the 
appearance of VO and ZrO features in the data, suggest that TP-AGB stars dominate these 
spectral features. However, a contribution from other evolved stars is also likely. Compar¬ 
ison with evolutionary population synthesis models shows that models based on empirical 
libraries that predict relatively strong NIR features provide a more accurate description of the 
data. However, none of the models tested accurately reproduces all of the features observed 
in the spectra. To do so, the models will need to not only improve the treatment of TP-AGB 
stars, but also include good quality spectra of red giant and E-AGB stars. The uninterrupted 
wavelength coverage, high S/N, and quantity of features we present here will provide a bench¬ 
mark for the next generation of models aiming to explain and predict the NIR properties of 
galaxies. 

Key words: galaxies: bulges — galaxies: active — galaxies: stellar content — stars: AGB 
and post-AGB 
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1 INTRODUCTION 


Studying the unresolved stellar content of galaxies generally in¬ 
volves disentangling the various components contributing to the 
spectral energy distribution (SED), fitting a combination of sim¬ 
ple stellar populations (SSPs) to derive information about age, 
metallicity, and star formation history. In the near-infrared (NIR, 
0.85-2.5/^m), the thermally pulsing asymptotic giant branch (TP- 
AGB) phase - the last stage of the evolution of intermediate mass 
stars (M<6M0) - is a particularly important component of the 
SSP models. These stars may be able to domin ate the emission 
of stell ar populations with ages ^ 0.2 - 2 Gyr Jm arasto ^| 2ooa 
see also lMouhcine & Lanconll2002h and lPottori et al. 1 2005h l be¬ 
ing responsible for roughly half of the luminosity in the K-band 
dSalaris et al.ll2014ll . 

Unfortunately, a correct treatment of the TP-AGB phase is 
difficult to achieve, since the phy sics of the evo l ution of this 
stellar phase is still poorly knovyn lMaras torj 200^Marigoet^^ 


2008l:IConroy & Gunrj20ia:IConrov et alj2013l:lKriek et al.l2010l: 


Zibetti et"^ 20131 : Noel et al] 2013h . The complex processes oc¬ 
curring during this phase (mass-loss, changing opacities, dredge-up 
events, etc.) are difficult to accu rately model, and m ust be anchored 
with empirical calibrations (e.g. lMarigo et al.ll2008l) . The contribu¬ 
tion of a population of TP-AGB stars to the light of a galaxy also 
depends on the duration of this evolutionary stage, which is not yet 
well constrained. In addition, the features of an observed spectrum 
of a TP-AGB star depend strongly on its age and metallicity (e.g. 
lUancon et alj|200i]) . and we do not have NIR empirical spectra of 
stars covering the full range of these parameters expected in galax- 


These difficulties have led to the production of models 
making very different predictions for the NIR spectral region. 
In contrast to Evolutionary Population Synt hesis (EPS) models 
based on isochrone synthesis, such as thos e of iB ruzual & Charlotl 
12003 , hereafter BC03), the models of iMarastonI j2005l) and 


iMaraston & StrdmbackI j201lL hereafter M05 and Mil), based on 
the fuel consumption theory and using empirical spectra of C- 
and 0-Rich stars for the TP-AGB phase, predict the presence of 
strong NIR molecular features, such as TiO (0.843 fj,m, 0.886 /rm), 
CN (1.1 /rm, 1.4 /rm), C 2 (1.17 ^m, 1.76 ^m), H 2 O (1.4 /rm, 
1.9 ^m) and CO (1.6 /rm, 2.3 /rm) at the ages where these stars 
are abundant. When applied to galaxy spectra, the M05 models 
also predict much hig her NIR luminosities than the BC03 ones 
jMaraston et al.l l2006h . The models are calibrated using Galac¬ 
tic and Magellanic Cloud globular clusters, and have been tested 


results iMarastonet^ ll2006l:lRiffeletaI 


2 OI 0 I: fciek et al.l 120101: iLvubenova et al 


2 OI 2 I: IZibetti et alj|2013l: lMarigoll2014l) . 


201 llJ: 

Lvubenova et al. 

I 2 OI 2 I: 

Melbourne et al. 


Given the challenges involved in constructing SSPs appropri¬ 
ate to modeling galaxy SEDs based solely on theoretical stellar 
spectra and/or limited empirical stellar spectra, searches for these 
bands in galaxies would be a valuable step forward in guiding the 
development and improvement of the models. Some studies have 
been carried out in this area. The 1.1 fj,m CN band , for exam¬ 
ple, has been detected in AGN an d starburst galaxies jRiffel et al.l 
l2007l : l^mos Alme ida et al.l 200^. alon g with the 1.4 fim CN band 
in a single object 1 Martins et al.ll20l'3^ . The CO bandheads near 
2.3 /rm are well known, and the first detection of ZrO features at 
0.8 - 1.0 Atnt was rec ently reported in a sample of star burst nuclei 
jMartins et alj|2013^ . Conversely, IZibetti et al.l (1201 3h do not de¬ 


tect the TP-AGB spectral features predicted by M05 in their spectra 
of post-starburst galaxies at z~0.2. 

Here we investigate the origin of the stellar molecular absorp¬ 
tion features in NIR spectra of 12 nearby galaxies, and use the spec¬ 
tra to test two current EPS models. The detected bands include CN 
(1.1 /rm, 1.4 fim), CO (>^1.6 fim, ~2.3 /rm), H 2 O (~1.4 /rm, ~1.9 
/^m), VO (1.048 ^tm) and TiO (0.843 /rm, 0.886 /rm). We also de¬ 
tect a feature at ~0.93 fim which is probably due to a blend of CN, 
TiO, and ZrO bands. A tentative detection of a C 2 (1.76 fim) band- 
head in the mean spectrum of the whole sample is also reported. 
While some of these bands have been discussed in previous work, 
the novelty of this study is the simultaneous detection and analysis 
of many features over a wide range of wavelengths. 


2 SAMPLE AND DATA 


2.1 The Sample 

The galaxies in this paper ( Table [p are a subse t of those observed 
for the Palomar XD project dMason et alj|201^ . That program ac¬ 
quired high-quality nuclear NIR spectra of galax ies in the Palomar 
survey of nearby galaxies dHo et al .11 19951 [199711 . using the Gem¬ 
ini Ne ar-Infrared Spectr ograph (GNIRS) on the Gemini North tele¬ 
scope dElias et al.ll2006ll . The GNIRS sample contains AGNs cov¬ 
ering a wide range of luminosity and accretion rate. Of these, a 
subsample was selected for this study on the basis of high signal-to- 
noise ratio (S/N) data, strong stellar continuum, weak AGN emis¬ 
sion, and excellent cancellation of the strong telluric absorption co¬ 
incident with the location of the CN and H 2 O bands around 1.4 
fim. Wh ile some molecular bands are evident in the remaining 38 
objects dMason et al.ll201^ . we restrict this discussion to the 12 
galaxies with the highest-quality data. 

It is important to realize that these are not necessarily the 
galaxies where features due to TP-AGB stars would be expected 
to be the strongest. That is, no selection due to stellar population 
properties was made. The galaxies are simply those in which the 
data permit a close examination of the various bands that have 
been suggested to be signatures of TP-AGB stars. As shown in 
Table [T] with the exception of NGC 2832 and NGC 2768, all 
of the galaxies are spirals. Optical stellar population studies have 
been carried out for several of the galaxies, and in most cases old 
ages are suggested (Table [TJ. Possible exceptions are NGC 2655, 
NGC 2768, and NGC 5005. The Hff line index in NGC2655 sug¬ 
gests a mean age of -^2 G\Tj|SiPchen^ 2006h. while studies dis- 
agree about whet her old dSil ch enkd 200^ ^ang et al.l l20()l) or 
intermediate-age dSerra et al.l200^: ICrocker et al.l20 0^ stars d om- 
inate in NGC 2768. In NGC 5005. ICid Fernandes et al.ld2004ah find 
that an intermediate-age stellar population contributes about 45% 
of the nuclear optical light. 

Inter-comparing stellar ages obtained from different spectral 
regions, apertures, and line indices/stellar population models can 
be highly problematic. Nonetheless we use this information in 0 
in which we draw some tentative conclusions about the type of stars 
causing the molecular absorption features observed in this galaxy 
sample. 


2.2 The data 

The spectra were obtained using GNIRS’ cross-dispersed mode, 
which provides simultaneous spectral coverage from ~0.8-2.5 fim. 
The 0 . 3 ” slit was used, providing R~1200, and was orientated at 
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the mean parallactic angle at the time of the observations. The ob¬ 
servations were executed in queue mode (Programs GN-2012B- 
Q-80, GN-2013A-Q-16) and the observing condition criteria al¬ 
lowed for observations to be taken with thin cirrus and seeing < 
1". The data were r educed using standard procedures, described in 
iMason et alj ( l2015h . 

Some stellar features of interest are located in regions of poor 
atmospheric transmission, in which the telluric cancellation pro¬ 
cess (dividing by a slightly shifted and scaled spectrum of an A 
star observed near in time and airmass) may leave artefacts in the 
spectrum. For instance, the 0.93 /rm feature lies close to telluric 
F120 absorption at ~0.95 /im. The authenticity of this feature was 
verified in two ways. First, the shifting and scaling applied to the 
standard star were adjusted, while verifying that the strength and 
profile of the feature were not significantly affected. Second, a stan¬ 
dard star spectrum was reduced using the same procedure applied 
to the galaxies, showing that no similar feature was produced in 
the resulting spectrum (Fig.[T). A similar procedure was applied to 
the H 2 O absorption regions, indicating that the broad F 120 bands 
apparent in some objects are not artefacts of the data reduction but 
true stellar features. 


3 RESULTS 

3.1 The Spectral Features 

The full set of spectra is shown in Figure[T] while Figures|2-|4]show 
close-ups of various regions of interest. As can be seen from these 
figures (see also Fig. [^, the CN (1.1 pim, 1.4 /rm) and CO (1.6 /rm, 
2.3 /rm) bands are clearly detected in the majority of so urces. While 
the 1. 1 jim CN band is known in a number of galaxies jRiffel et al.l 
l2007h . the CN band at 1.4 ixm h as to date only been detected in 
NGC 4102 dMartins et al.ll2()13ah . The CN 1.4 ^m band falls in a 
region of poor atmospheric transmission and is therefore difficult to 
detect in low-redshift extragalactic sources. Thanks to the dry con¬ 
ditions under which these data were obtained, however, we detect 
the 1.4 /rm CN band in almost all the galaxies of the present sample. 
According to M05, these bands are strong in TP-AGB stars, whose 
contribution to the integrated light peaks at intermediate ages. 

The FI 2 O 1.4 /rm band is also detected in almost all sources, 
while the 1.8 ^m band is visible only in some (NGC 2655, 
NGC 2768, NGC 3147, NGC 4548, NGC 5850). C 2 lines are weak 
or absent in the individual galaxy spectra but the 1.76 pm band 
may show up in the mean spectrum (see Fig. O- These features 
are present in the IRTF a tlas of cool stars ( s ee Fig. 7 to 34 of 
iRavner et alj2009L see also lCushing et alj l l2005h ). and are also pre¬ 
dicted by M05 (their Figures 14 and 15). 

A comparison of the GNIRS spectra with the IRTF library 
shows that, besides the bands predicted by the M05 models, we 
also detect a VO feature at 1.048 pm. A comparison of our Fig.[T] 
with Fig. 8 of iRavner et all ( l2009l) shows that the VO band at 1.048 
pm, detected in many of the galaxies in our sample, most likely 
originate in late M III stars in the TP-AGB phase. To our knowl¬ 
edge, this is the first time that these bands have been reported in 
galaxies. A broad feature ranging between ^0.93 - 0.95 pm is also 
detected, but the identity of the carrier is uncertain. The feature may 
be due to a combination of CN (red system, Au = v” —v' = —1), 
TiO (e system, Av — — 1), and ZrO, which are all strong in the 
spectra of cool giant stars jRavner et alj200^ . The 0.886 pm band 
due to the S system of TiO (An = 0) is detected in some of the 
spectra (NGC 2832, NGC 4594, and perhaps others). Depending 


Table 2. Line limits and continuum bandpasses. 


ID 

line limits 
(/ 2 m) 

continuum bandpass 
((im) 

CN/TiO/ZrO 

0.9200-0.9500 

0.8760-0.8800,0.9570-0.9650 

VO 

1.0470-1.0650 

1.0350-10390, 1.0660-1.0740 

CNll (q) 

1.0900-1.1300 

1.0660-10740,1.1310-1.1360 

CN14a 

1.4050-1.4300 

1.3950-13990, 1.4340-1.4390, 1.4750-1.4840 

CN14b 

1.4500-1.4715 

1.3950-13990, 1.4340-1.4390, 1.4750-1.4840 

CO 16a 

1.5735-1.5810 

1.5110-15170, 1.5390-1.5410, 1.6270-1.6310, 1.6570-1.6580 

C016b(cK) 

1.6175-1.6285 

1.5110-15170, 1.5390-1.5410, 1.6270-1.6310, 1.6570-1.6580 

C022a(Q) 

2.2910-2.3200 

2.2690-22790, 2.3950-2.3999 

C022b(a) 

2.3200-2.3420 

2.2690-22790, 2.3950-2.3999 

C022c(q) 

2.3420-2.3695 

2.2690-22790, 2.3950-2.3999 


Table Notes: a based on iRiffel et al.l l2009h and iRiffel et alj 1201 Ibl) . 


on the spectral coverage, the 0.843 pm An = 0 bandheads of the e 
system of TiO are also detected in those objects, blended with the 
well-known Calcium triplet lines. 

We have measured the equivalent widths (EW) of the strongest 
molecular bands, including the 0.93 pm CN/TiO/ZrO and 1.048 
pm VO bands that are particularly relevant to TP-AGB stars. The 
measurement s were performed using an updated version of the 
PACCE code jRiffel & Valell20l7l) . In this code, the continuum un¬ 
certainty is taken as the root mean squere of the difference be¬ 
tween the linear fit to the pseudo continuum and the observed one. 
Then a Gaussian distribution of errors, using the above difference 
as 1-sigma standard deviation, is assumed to simulate new pseudo- 
continua. Finally, EW uncertainties are assumed to be the standard 
deviation of 500 measurements of EW using the simulated pseudo 
continua as continuum points for the linear fit. The line definitions 
used are listed in Tablel^and the measured values are in Tablej^ 

As the spectra of the 12 galaxies have rather similar contin¬ 
uum and emission line properties (see Fig.[T}, we averaged them to 
construct a spectrum with very high S/N, allowing further verifica¬ 
tion that the detected features are real. The spectra were normalized 
to unity at 1.223 pm before computing the mean value of the fluxes, 
pixel-by-pixel. The resulting spectrum (Fig.[3 supports the conclu¬ 
sion that the features observed in the individual galaxies are real, 
since they clearly appear in the mean spectrum. As the galaxies are 
at different redshifts, when they are shifted to the rest frame to cre¬ 
ate the mean spectrum, any systematic effects are shifted to differ¬ 
ent positions and will therefore tend to cancel out or be weakened 
in the mean spectrum. Additionally, (i) reducing a standard star us¬ 
ing the same procedures results in a very clean spectrum with few 
weak telluric line residuals (^, and (ii) we can fit the observed 
mean absorption spectrum very well using the IRTF stellar library 
(S|4j, thus indicating that the features are real and due to stars. 


3.2 Stellar population synthesis 


Bands from carbon-containing molecules are present in the galaxy 
spectra in Fig. [T] along with VO and possibly a contribution from 
ZrO in the band at 0.93 pm. The presence of these molecules sug¬ 
gests a significant contribution from TP-AGB stars to the NIR spec¬ 
tra of the galaxies. To further investigate (1) whether these molecu¬ 
lar bands really are due to TP-AGB stars, and (2) whether the avail¬ 
able stellar population synthesis models are able to reproduce them, 
we carry out a stellar population synthesis o n the mean spectrum 
using the STA R LIGHT spectral sy nthesis code jCid Fernandes et ^ 
l2004aLl2005all5:lAs ari et alj200 ^. The procedure s followed are de¬ 
scribed in Riffel et al. ( 2009h and lPametto et al.l j2014l) . including. 


for example, the methods used to handle extinction, emission lines 
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Figure 1. Rest frame spectra of the galaxy sample. Shaded regions indicate broad, stellar H 2 O molecular bands, while dotted lines denote the CO bands. 
Other molecular bands are labeled at the bottom of the figure. A Mauna Kea atmospheric transmission spectrum for 1 mm H 2 O and 1 airmass is also shown 
(gray), along with the spectium of an AO stai* (blue) reduced in the same manner as the galaxy data. Because of the small redshifts of the galaxies, the telluric 
features/residuals in the AO star and atmospheric spectra do not correspond precisely to those in the galaxy data. 
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Table 1. Galaxy Properties 


Galaxy 

Distance^* 

Mpc 

Morphology" 

AGN" 

NGC 2655 

24.4 

Sa 

S2 

NGC 2768 

23.7 

S0/E6 

L2 

NGC 2832 

91.6 

E3 

L2 

NGC 3147 

40.9 

S(s)b I-II 

S2 

NGC 3718 

17.0 

Sa 

LI.9 

NGC 4548 

16.8 

SB(sr)b I-II 

L2 

NGC 4565 

9.7 

Sb 

S1.9 

NGC 4594 

20.0 

Sab 

L2 

NGC 4750 

26.1 

S(r) 

LI.9 

NGC 5005 

21.3 

S(s)b II 

LI.9 

NGC 5371 

37.8 

SB(sr)b I 

L2 

NGC 5850 

28.5 

SB(rs)b I-II 

L2 


log L(Ha)°‘ 

Age (ref)*’ 

Aperture" 

SNR'! 

erg s“^ 

39.55 

1(1) 

pc 

27x82 

118 

39.01 

1(1,2), 0(3,4) 

27x80 

121 

38.91 

O(2,5) 

135x404 

160 

39.47 

0(6) 

54x163 

132 

38.46 

- 

19x58 

123 

38.46 

0(7) 

9x28 

no 

37.97 

- 

24x73 

100 

39.70 

0(8) 

21x63 

168 

39.04 

- 

31x94 

119 

39.47 

1/0(9) 

19x56 

119 

39.24 

- 

50x149 

100 

38.66 

0(10) 

50x149 

100 


Table Notes: (a) 


Ho et alj jl997b and ref ere nces therein, (b) S tell ar population age. O = ol d, age >5 Gyr. I = in ter mediate-age, age ^1-2 Gyr. Re ferences: 1 


ISil’cheni^ flOOf 

<2005t) . 8-1^ 

X 0''.9. (d) Determined in the range 2.076-2.096 /rm. 


2-IZhang et all feOOSh 3-ls iH200a. 4-1 I ^C^cker et al. 

Sanchez-Blazquez et al. i2006l) . 9- ICid Femand^et al 


10 - 


llQOSh. 5-lLoubser et al, 


de Lorenzo-Caceres et alj 1201 


6- lBruzual & Chariot I l l2003h . 7- ISai'zi et alJ 
(c) Slit width X extraction aperture of 0''.3 


Table 3. Equivalent widths of the absorption bands in the galaxy sample. 


Galaxy 

CN/TiO/ZrO 

VO 

CNll 

CN14a 

CN14b 

CO 16a 

CO 16b 

C023a 

C023b 

C023c 

NGC 2655 

20.20±0.50 

3.66±0.37 

12.06±0.52 

8.31±0.50 

4.71 ±0.34 

4.12±0.13 

6.84±0.19 

23.31±0.24 

22.60±0.30 

29.83±0.42 

NGC 2768 

17.23±0.53 

3.75±0.24 

1L14±0.46 

5.09±0.58 

5.88±0.37 

3.23±0.07 

6.45±0.11 

23.97±0.27 

22.61±0.22 

30.42±0.33 

NGC 2832 

16.16±0.49 

2.17±0.35 

1L24±0.51 

8.88±0.34 

5.00±0.27 

4.30±0.10 

4.85±0.27 

21.59±0.43 

2L24±0.43 

30.63±0.68 

NGC 3147 

23.80±0.41 

3.23±0.46 

9.45±0.51 

6.68±0.42 

3.62±0.38 

2.39±0.08 

5.87±0.16 

18.92±0.21 

17.44±0.29 

24.18±0.40 

NGC 3718 

18.34±0.58 

4.92±0.32 

9.07±0.61 

8.22±0.49 

4.92±0.39 

3.37±0.09 

6.28±0.19 

22.36±0.33 

20.58±0.26 

27.90±0.37 

NGC 4548 

22.58±0.75 

3.94±0.45 

10.49±1.15 

6.33±0.70 

4.74±0.47 

3.33±0.10 

6.90±0.17 

20.33±0.25 

18.77 ±0.25 

25.97±0.36 

NGC 4565 

16.61±0.44 

3.43±0.24 

9.01±0.49 

7.63±0.49 

4.16±0.31 

3.64±0.08 

6.66±0.17 

21.98±0.25 

20.94±0.24 

28.45±0.42 

NGC 4594 

20.07±0.40 

4.81 ±0.24 

8.07 ±0.40 

10.43±L02 

3.74±0.37 

3.41 ±0.07 

3.68±0.16 

23.80±0.16 

22.13±0.16 

30.86±0.32 

NGC 4750 

16.44±0.61 

3.98±0.40 

8.32±0.56 

5.37±0.55 

4.43±0.33 

L96±0.08 

6.41±0.15 

20.86±0.28 

20.61±0.30 

26.45±0.48 

NGC 5005 

11.86±0.41 

L70±0.32 

10.52±0.53 

7.93±0.23 

3.39±0.27 

4.28±0.08 

6.15±0.13 

22.80±0.22 

23.16±0.26 

30.06±0.48 

NGC 5371 

16.80±0.48 

L75±0.34 

7.26±0.67 

4.84±0.38 

3.81±0.35 

L94±0.12 

6.41±0.21 

22.65±0.25 

20.62±0.28 

28.82±0.44 

NGC 5850 

18.45±0.51 

3.11±0.39 

10.06±0.66 

3.85±0.44 

3.87±0.31 

3.07 ±0.09 

6.72±0.14 

22.58±0.30 

2L94±0.26 

29.15±0.34 


and differences in spectral resolution. Briefly, STARLIGHT fits an 
observed spectum 0\ with a combination, in different proportions, 
of SSPs, solving the equation: 


M\ = Mao 


■ iV* 

J=l 




( 1 ) 


where M\ is a model spectrum, r\ is the reddened spectrum 
of the jth SSP normalised at Ao; r\ = is the red¬ 

dening term; Maq is the theoretical flux at the normalisation wave¬ 
length; X is the population vector, and G(n,,( t,) is the Gaussian 
distrihution used to model the line-of-sight stellar motions, which 
is centred at velocity n* with dispersion cr*. The final fit is carried 
out searching for the minimum of the equation: 


A 


( 2 ) 


tio n of STARLIGHT see i t s manu afi For applications i n the optical 
see ICid Fernandes et alj (l2005all . and in the NIR see iRiffel et Tl 

(l20m 

We use 4 different approaches to define the base set of spectra 
used by the code, as follows: i) stars approach, a base composed 
of all the 210 dereddened stars (spectral types F-S/C, F being the 
hottest available) in the IRTF spectral library jRavner et al.ll200^ : 
ICushing et al.ll20(i^ and ii) stars no TP-AGB approach: the same 
as in (i) but removing all the TP-AGB stars from the base, hi) Mil 
approach: EPS models of Mil (an example of models including 
empirical spectra of C- and O-Rich TP-AGB stars), which consist 
of theoretical SSPs covering 12 ages {t =0.01,0.03, 0.05, 0.1, 0.3, 
0.5, 0.7, 1, 2, 5, 9 and 13 Gyr) at solar metallicity; iv) BC03 ap¬ 
proach: SSPs from BC03 (an example of models calibrated with 
broad band photometry) with the same ages and metallicity as 
above. 


where emission lines and spurious features are masked out by fix¬ 
ing w\=0. The quality of the fit is accessed by X/Sed wich is the 
given by equation divided by the number of points used in 
the fit and by adev = |Oa — Ma|/Oa, which is the percent¬ 
age mean deviation over all fitted pixels. For a detailed descrip- ^ http://astro.ufsc.br/starlight/ 
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Figure 2. Close-up of the 0.8 — 1.6 /xm region. Top: Mil, M05 and BC03 models t l|3.2t . Other panels: observed spectrum (thin red line) and degraded to 
M05/BC03 spectral resolution (thick black line). 


4 DISCUSSION 

4.1 Origin of the features 

To link the molecular absorption features to the kind of star in 
which they may arise, we first revisit some relevant aspects of the 
later stages of stellar evolution. 

Once the helium in the core of a star is exhausted. He burn¬ 
ing begins in a shell surrounding the core. This marks the arrival 
of the star on the asymptotic giant branch (AGB). The star starts 
to expand, becoming cool and luminous, and develops a growing 
convective outer region. At this stage the second dredge-up occurs, 
mixing the end points of helium burning, mainly C and O, into 
the star’s atmosphere. At the end of this phase the s process be¬ 
comes the dominant me chanism of nuclear fusion , producing ele¬ 
ments such as Zr and V dHabing & 01ofss^l2004l) . 

The star then starts to undergo thermal pulses, leaving the 
E-AGB stage and entering the TP-AGB phase. The third dredge- 


up takes place, and processed material from the star’s interior is 
transported into its atmosphere. At this point the spectrum of the 
star will present evidence of enhanced C and O abundance, along 
with features due to specie s containing the newly-formed Zr and V 
dHabing & 01ofssonll2004) . For example, the NIR CN absorption 
band is particularly enhanced in carbon stars, where there is resid¬ 
ual carbon that is not bound with oxygen in CO molecules. As N 
type carbon stars are prod uced after the third convec tive dredge-up 
along the TP-AGB phase dHabing & 01ofssoiill2004l) . the detection 
of deep CN features may indicate stellar populati ons rich in this 
kind of star dUancon et alJIZOoH : iRiffel et al.ll2007h . Similarly, VO 
and ZrO are produced with the residual oxygen not used in the for¬ 
mation of CO and TiO. Since V and Zr are only produced by the s 
process, these features are normally only detected in TP-AGB stars. 

Below, we relate the spectral features in the galaxies to those 
observed in various kinds of star. It is therefore necessary to un¬ 
derstand the stellar characteristics (spectral types and other traits) 
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Figure 3. Same as figure[2] 


that identify a star as being in the T P-AGB phase. Most AGB stars 
have spectral classes M, S, and C jHabing & 01ofssoi]l2004ll . Of 
the two main types of carbon star, C-N stars are in the TP-AGB 
phase while C-R stars are not. All the molecular bands in C stars 
are from carbon-containing molecules. The spectra of early M-type 
stars contain TiO bands, while VO becomes apparent at later spec¬ 
tral types. S stars, however, are characterized by distinctive ZrO 
absorption. As noted above, s-process elements such as Zr and V 
are normally only present in AGB stars in the thermally-pulsing 
stage. 

TP-AGB stars can also be distinguished by the variabil- 
ity caused by their long-period pulsations. The stellar library of 
iRavner et alJ ^20091) includes several IR spectra of TP-AGB stars, 
identified by their characteristic variability types (irregular, semi¬ 
regular, and Mira), and their data illustrate the above trends. In¬ 
specting their series of dwarf, giant and supergiant spectra (their 
Figures 7, 8, and 9) shows that the only molecular features detected 
in the dwarf stars are TiO, CO, HgO and FeH. CN bands are present 


in the giants and supergiants, but are most prominent in the late- 
type supergiants with TP-AGB-like variability characteristics (the 
K1 la-Iab, K5 Ib, M3 lab-Ia, and M5 Ib-II stars in their Figure 9). 
ZrO also becomes prominent in the coolest of the TP-AGB super¬ 
giants in that figure. In the latest-type of the TP-AGB giants (M4 
III, M6 III, M7 III, and M8-9 III in Figure 8 of Rayner et al. 2009), 
absorption from the VO molecule becomes apparent. 

The results of the stellar population synthesis are shown in 
Figure]^ When using the stars approach the spectrum can be well 
reproduced, and a significant fraction (~ 35%) of the light comes 
from TP-AGB stars (the stars marked with * in the lower left panel 
of the figure). This is illustrated by Figured which shows that 
these stars account for the molecular features particularly at 
the short wavelength end of the spectrum. However, when the 
TP-AGB stars are removed, the code replaces them with a variety 
of E-AGB and red giant stars, notably a large contribution from a 
K5Ib star and a significant contribution from a C-R type carbon star. 
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Figure 4. Same as figure[2] 


The quality of the fit is similar overall, with just a small increase in 
the and Aderll values. 

Taken at face value, this would suggest a degeneracy between 
R-type carbon rich stars summed with E-AGB stars, and the TP- 
AGB stars. However, there are two important details that have to 
be taken into account. First, we are trying to fit a galaxy’s stellar 
population with the observed spectra of stars from a library which 
is neither complete nor absolutely flux-calibrated, and most im¬ 
portantly, not constrained by the rules of stellar evolution, which 
would require the use of an isochrone and an assumption about the 
initial mass function.This would be the correct approach to cre¬ 
ate an ssilE where the contribution of each stellar spectral type 
is weighted by the number and luminosity of the stars. This is a 
large task in its own right and is well beyond the scope of this pa¬ 
per. This means that although we can reproduce the mean spectrum 
with both approaches, the results need to be taken with caution. 
AGB and regular giant stars are short-lived, and their true contri- 


^ The percentage mean | I /0;\ deviation over all fitted pixels. 

® Which are the coiTect inputs for the base of elements in stellar population 
fitting using the STARLIGHT code. 


bution to a stellar population spectrum of a given age cannot be 
determined by the approaches we described here. 

With that said, TP-AGB stars are very luminous. Despite the 
fact that they exist only in population components with ages ~0.2 
- 2 Gyr (M05), they have been show n to account for almos t 20% 
of the 1.6 ^m flux of nearby galaxies l lMelboume et alj2012h . with 
much of the remainder likely coming from relatively featureless 
main sequence stars. In the K-ba nd they may be re sponsible for 
roughly half the total luminosity l lSalaris et alj|2014h . Such lumi¬ 
nous stars, with strong absorption features, may be expected to 
leave an imprint on the NIR spectrum. Also, there is a subtle dif¬ 
ference between the two fits in the absorption band at 0.93 ^m (see 
Fig.0. As noted above, this band probably contains contributions 
from CN, TiO, and ZrO. The fit in this region is improved when 
TP-AGB stars are included, suggesting that some of the absorp¬ 
tion indeed comes from the Zr brought to the surface by the third 
dredge-up. Combined with the fact that VO is present in the galax¬ 
ies’ spectra (see Fig.[^ these considerations suggest that TP-AGB 
stars pl^ a significant role in the features observed in our integrated 

spectrqS 

^ We note that the VO feature near 1.05 fim is not well reproduced by 
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Wavelength (/xm) Wavelength (/xm) Wavelength (/xm) 



Figure 5. Top: Mean spectrum (black), along with EPS fits using (lowest to highest) the IRTF stars, stars no TP-AGB, Mil and BC03 models. Middle: the 
difference between the mean spectrum and its fits, following the same order a s in the top. Bottom: s ynthesis results from these three methods. Xj denotes 
light fractions and denotes mass fractions. Stars spectral types are taken from ICushing et aO j2005h : lRavner et alJi2009h and * denotes TP-AGB stars. The 
goodness of the fits and Av are on the labels. For display purposes we removed from the histograms stars in the base set that do not contribute to the fits. 


If TP-AGB stars dominate the NIR spectral features, we might 
expect to observe stronger features in intermediate-age galaxies 
that should host larger populations of these stars. This is not ob¬ 
viously the case in this sample. The spectra in Figure|T]are remark¬ 
ably similar, and there is no clear relationship between the EW val¬ 
ues in Table [3 and whether a galaxy is classified as intermediate- 
age or old in Table [T] Of the galaxies in this sample, NGC 5005 
has probably the most secure ev idence for a large intermed iate- 
age (-^45%) stellar population jCid Fernandes et all l2004bh . Its 
CN/TiO/ZrO and VO bands, however, are the weakest that we 
observe (Table [3. Similarly, iMason et alj l l2015h examined three 
galaxies with intermediate-age populations and found them to have 
weaker molecular absorption bands than older galaxies. 

There are several possible reasons for this apparent contradic¬ 
tion. First, stellar populations derived from optical and NIR spec¬ 
troscopy cannot necessarily be directly compared. Dust can obscure 
young populations in the optical, and different slit size!0 can probe 


either the models or the stellar library. The absence of the feature in the 
model spectra is not surprising, as they do not include the VO molecule. Its 
weakness in the empirical fit may be a result of incompleteness of the stellar 
library or the fact we are not constraining it by the rules of stellar evolution. 
A definite conclusion will require more detailed analysis, and is left for a 
future paper (R. Riffel et al in preparation) 

® That might be the case of NGC5005, since the slit used by 


different regions of s patially varied populations l lRiffel et aljTOlld: 
iMartins et alj|2013bl). Second, the intermediate-age galaxies stud¬ 
ied here and in Mason et al. I l20l^ likely span only limited range 
of metallicity. As the strength of TP-AGB spectral features depends 
greatly on their age and metallicity (M05), weak NIR features like 
those observed in NGC 5005 may not be a universal characteristic 
of intermediate-age galaxies. Third, intermediate-age galaxies may 
also contain young stellar populations, whose relatively featureless 
spectra may dilute the stronger bands of the TP-AGB stars. 

Finally, as shown by our model fits, red giants, C-R and E- 
AGB stars can also produce most of the observed features some¬ 
times attributed to TP-AGB stars. These stars may be important in 
stellar populations that are both younger than the intermediate-age 
populations in which TP-AGB stars are believed to be abundant 
(if R C-rich and E-AGB stars are present in large fractions), or in 
older populations (which contain large numbers of red giants). This 
clearly demonstrates that to accurately model a stellar population, 
all kinds of evolved stars need to be taken into account. 

To date, the most targeted s earch for TP-AGB stars in galax¬ 
ies is that of IZibetti et alj l l2013h . They obtained NIR (H and K- 
band) spectra of 16 post-starburst galaxies, selected partly for the 
CN (1.41 /rm) and C 2 (1.77 ^m) bands to be redshifted out of 


ICid Fernandes et ^ l2004bh is three times wider than ours, thus potentially 
sampling additional stellar populations. 
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a(A) 


Figure 6. Fractional contribution of the stars used in the fit using the stars approach (Top). For display purposes we removed the stars contributing with less 
than 2% (to see the contribution of stars contributing with less than 2% see Fig.|^ from the plot. The contribution of the TP-AGB stars (marked with *) to the 
observed absorption features is clearly observed, a zoom to these features is shown at the bottom. For more details see Sec l3.2l 



Figure 7. Difference between the stars and stars no-TP-AGB approaches 
presented in Fig. [5] 


strong telluric absorption (z~0.2). Although the stellar ages (0.5 
-1.5 Gyr) of these galaxies correspond to the period of maximum 
contribution of TP-AGB stars to the integrated spectrum, the depth 


of absorption features from carbon molecules in these stars was 
found to be no more than a few per cent of the continuum. Ex¬ 
trapolating from our observations to predict the features expected 
in post-starburst galaxies is beyond the scope of this paper, but we 
note that the weak CN and C 2 absorptions reported here would not 
be detected at the resolution (^300 versus 1200) and S/N (-^30 
vs 110 in the K band) of the Izibetti et al.l l2013h data, and their 
spectral coverage does not include most of the other features that 
we report. It is therefore difficult to gauge whether our interpre¬ 
tation of the molecular features as dominated by TP-A GB stars is 
in conflict with the observations of IZibetti et alj ( l2013h . However, 
it is clear that the features in the post-starburst galaxies are not as 
strong as predicted by the M05 models for galaxies of those ages 
and metallicities. 

4.2 Performance of the models 

The last two approaches tested in 113.21 (Mil approach and BC03 
approach) are intended to show how well different SSP recipes are 
able to reproduce these features. Fitting the mean spectrum with the 
Mil models results in a complex star formation history, with simi¬ 
lar fractions of young (t < 1 Gyr), intermediate-age (t ~ 1 Gyr) and 
old (t ~ 10 Gyr) SSPs. With the BC03 models, though, STARLIGHT 
includes a young (~ 10^ yr) component and requires only a very 
small fraction of old SSPs. This is probably because the BC03 mod¬ 
els do not have deep NIR CN bands, causing the code to attempt 
to fit the 1.1 /rm CN band with Pa 7 . If such a young component 
were genuinely present, we would expect to require a significant 
fraction of F stars (the hottest available in the IRTF library) when 
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fitting the mean spectrum with the stellar library, which is not the 
case. Furthermore, if we remove the ^100 Myr components in the 
BC03 approach, the best fit is achieved with ~ 35% of the light due 
to the 1 Gyr and 65% due to the 2 Gyr components, but with a 
much poorer fit (x^ ~ 100). 

The high (~65%) light fraction of stars with ages 0.2-3 x 10® 
yr in the Mil model fits ostensibly supports our previous interpre¬ 
tation that the molecular bands likely contain a strong contribution 
from TP-AGB stars, which should dominate the NIR emission of 
a stellar population of this age. However, this result could be co¬ 
incidental; the Mil models include empirical spectra with good 
spectral resolution only for the TP-AGB stars, meaning only the in¬ 
termediate age stellar population in their models will present these 
features. As discussed above, most of the features can be repro¬ 
duced with a mix of E-AGB and R C-rich stars, which means older 
stellar populations probably also have these features at some level. 
To avoid misinterpreting NIR galaxy spectra, it is thus of great im¬ 
portance to include high quality empirical spectra of the other evo¬ 
lutionary phases in the EPS models. 

While the Mil models qualitatively reproduce several of the 
broad absorption bands in the mean spectrum (notably CO and 
CN), these bands are weaker or absent in the BC03 model fit. This 
is probably not related to the different spectral resolution of the 
models (R ~ 200, vs i? ~ 500 in Mil) because we find that 
similar results are obtained from the Mil and M05 models (the lat¬ 
ter with R 250 in the K-band). The differences are more likely 
related to uncertainties in the theoretical treatment of AGB and TP- 
AGB stellar phases and the fact that the M05/M11 models include 
empirical spectra of C and O-rich stars and therefore better account 
for these features. 

Since the relative band strengths depend on metallicity (M05), 
we also tried fitting the mean spectrum with BC03 and M05 mod¬ 
els including three non-solar metallicities (the Mil models include 
only solar metallicity for this wavelength range). Similar ages to 
those derived from the original analysis were obtained in all cases, 
and changing the metallicity did not significantly improve the fits 
to the features. This may suggest that the incompleteness in stel¬ 
lar types/phases of the empirical spectra on which the models are 
based is the first order problem in reproducing the NIR absorption 
features. At this stage in the development of the models, metallic¬ 
ity probably plays a secondary role in predicting the strength of the 
features 


5 FINAL REMARKS 

We have presented high-S/N, 0.85 - 2.5 ^m spectra of the nuclei of 
12 local galaxies. The spectra show indications of numerous weak 
molecular absorption bands that arise in the atmospheres of cool 
stars. While several types of evolved stars can exhibit most of the 
detected features, we argue that the high luminosity of TP-AGB 
stars, and the fact that the spectra contain VO and likely ZrO, that 
are only expected in stars undergoing thermal pulses, imply that 
these features are the long-debated signatures of TP-AGB stars. 
However, a combination of R-type carbon rich stars and E-AGB 
stars should not be ruled out as an important contributor to these 
bands. 

In testing two sets of stellar population synthesis models 
against these spectra, we found that models that include stars 
with relatively strong absorption bands provide a more accurate 
representation of the spectra. However, to avoid misinterpreting 
the data, future models must include high-quality spectra of C-R 


and E-AGB stars, along with those of TP-AGB stars. Testing the 
predictions of EPS models has to date been hampered by a lack 
of high-quality data againstwhich the models can be tested, and 
we hope that the high S/N and wavelength coverage of the spectra 
presented here will make them a useful aid to the development 
of such models. The reduced spectra presented here are available at: 
http://www.canfar.phys.uvic.ca/vosui/#/karun/xdgnirs_Dec2014/XD .final 


ACKNOWLEDGEMENTS 

We thank the anonymous referee for useful comments. The Brazil¬ 
ian authors thank CNPq and FAPERGS support. RR acknowl¬ 
edges M. Cushing, L. Girardi, C. Maraston, A. Romero and R. P. 
Schiavon for helpful discussions as well as CNPq and FAPERGS 
for financial support. LCH acknowledges support by the Chinese 
Academy of Science through grant No. XDB09030102 (Emer¬ 
gence of Cosmological Structures) from the Strategic Priority Re¬ 
search Program and by the National Natural Science Foundation of 
China through grant No. 11473002. CRA is supported by a Marie 
Curie Intra European Fellowship within the 7th European Commu¬ 
nity Framework Programme (PIEF-GA-2012-327934). Based on 
observations obtained at the Gemini Observatory, operated by the 
Association of Universities for Research in Astronomy, Inc., un¬ 
der a cooperative agreement with the NSF on behalf of the Gemini 
partnership: the National Science Foundation (US), the Science and 
Technology Facilities Council (UK), the National Research Coun¬ 
cil (Canada), CONICYT (Chile), the Australian Research Council 
(Australia), Ministerio da Ciencia e Tecnologia (Brazil) and Minis- 
terio de Ciencia, Tecnologia e Innovacion Productiva (Argentina). 


REFERENCES 

Alonso-Herrero, Almudena; Rieke, Marcia I.; Rieke, George H.; 

Shields, Joseph C., 2000, ApJ, 530, 688. 

Asari,N. V., Cid Fernandes, R., Stasihska, G., Torres-Papaqui, J. 
R, Mateus, A., Sodre, L., Schoenell, W., Gomes, J. M., 2007, 
MNRAS, 381, 263 

Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Cid Fernandes, R., Gu, Q., Melnick, J., Terlevich, E., Terlevich, 
R., Kunth, D., Rodrigues Lacerda, R., Joguet, B., 2004, MN¬ 
RAS, 355, 273 

Cid Eemandes, R.; Gonzalez Delgado, R. M.; Schmitt, H.; 
Storchi-Bergmann, T.; Martins, L. R; Rerez, E.; Heckman, T.; 
Leitherer, C.; Schaerer, D., 2004, ApJ,605,105. 

Cid Eemandes, R., Mateus, A., Sodre, Laerte, Stasinska, G., 
Gomes, J. M., 2005a, MNRAS, 358, 363 
Cid Fernandes, R., Gonzalez Delgado, R. M., Storchi-Bergmann, 
T., Martins, L. Fires & Schmitt, H., 2005b, MNRAS, 356, 270 
Chariot, Stephane; Worthey, Guy; Bressan, Alessandro, 1996, 
ApJ, 457, 625. 

Chisari, Nora E.& Kelson, Daniel D., 2012, ApJ, 753, 94. 
Crocker, A. F, Bureau, M., Young, L. M., Combes, F. 2008, MN¬ 
RAS, 386, 1811 

Conroy, Charlie; Gunn, James E., 2010, ApJ, 712, 833. 

Conroy C., 2013, ARA&A, 51, 393. 

Cushing, M.C. Rayner, J.T., & Vacca, W.D., ApJ, 2005, 623, 1115 
Dametto, N. Z.; Riffel, R.; Fastoriza, M. G.; Rodn'guez-Ardila, 
A.; Hernandez-Jimenez, J. A.; Carvalho, E. A., 2014, MNRAS, 
443, 1754. 


© RAS, MNRAS 000,[T]-?? 


12 Rijfel et al. 


de Lorenzo-Caceres, A., Falcon-Barroso, J., Vazdekis, A. 2013, 
MNRAS, 431,2397 

Dottori, H.; Diaz, R. J.; Carranza, G.; Llpari, S.; Santos, J., Jr. 
Elias, Jonathan H.; Joyce, Richard R.; Liang, Ming; Muller, 
Gary R; Hileman, Edward A.; George, James R., 2006, SPIE, 
6269E,4CE. 

Gonzalez Delgado, Rosa M.; Cid Fernandes, Roberto; Pfez, En¬ 
rique; Martins, Lucimara R; Storchi-Bergmann, Thaisa; Schmitt, 
Henrique; Heckman, Timothy; Leitherer, Claus, 2004, ApJ, 605, 
127. 

Girardi, Leo; Marigo, Paola; Bressan, Alessandro; Rosenfield, 
Philip, 2013, ApJ, 777, 142 

Habing, H. J & Olofsson, H., 2004, Springer, ISBN: 978-1-4419- 
1843-7 

Ho, L. C.; Filippenko, A. V.; Sargent, W. L., 1995, ApJS, 98, 477. 
Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997, APJS, 
112, 31 

Ho, L. C.,2008, ARA&A, 46, 475. 

Hunt, L. K.; Thuan, T. X.; Izotov, Y. L, 2003, ApJ, 588, 281. 
Kelson, Daniel D.; Holden, Bradford R, 2010, ApJ, 713, 28. 
Kriek et al., 2010, ApJ, 722L, 64. 

Kuntschner et al., 2010, MNRAS, 408, 97. 

Lan^on, A. & Wood, P. R.Goldader, 2000, A&AS, 146, 217. 
Lanjon, A., Goldader, J. D., Leitherer, C., & Gonzalez Delgado, 
R. M., 2001, ApJ, 552, 150 

A. Lanfon, P. H. Hauschildt, D. Ladjal & M. Mouhcine, 2007, 
A&A, 468, 205. 

Loubser, S. I., Sanchez-Blazquez, R, Sansom, A. E., Soechting, 1. 

K. 2009, MNRAS, 398, 133 

Lyubenova M., Kuntschner H., Rejkuba M., Silva D. R., Kissler- 
Ratig M., Tacconi-Garman L. E., Larsen S. S., 2010, A&A, 510, 
A19 

Lyubenova M., Kuntschner H., Rejkuba M., Silva D. R., Kissler- 
Ratig M., Tacconi-Garman L. E., 2012, A&A, 543, A75 
Mason, R. E. et al., 2015, ApJS, 217, 13. 

Maraston, C., 2005, MNRAS, 362, 799. 

Maraston, C.; Daddi, E.; Renzini, A.; Cimatti, A.; Dickinson, M.; 

Rapovich, C.; Rasquali, A.; Rirzkal, N., 2006, ApJ, 652, 85. 
Maraston, C. & Strbmback, G., 2011, MNRAS, 418, 2785. 
Marigo R, Girardi L., Bressan A., Groenewegen M. A. T, Silva 

L. , Granato G. L., 2008, A&ASS, 482, 883. 

Marigo R, 2014, arXiv:141L3126 

Martins, Lucimara R; Riffel, Rogerio; Rodriguez-Ardila, Alberto; 

Gruenwald, Ruth; de Souza, Ronaldo, 2010, MNRAS, 406, 2185 
Martins, Lucimara R; Rodriguez-Ardila, Alberto; Diniz, Suzi; 
Gruenwald, Ruth; de Souza, Ronaldo, 2013, MNRAS, 431, 
1823. 

Martins, Lucimara R; Rodriguez-Ardila, Alberto; Diniz, Suzi; 

Riffel, Rogerio; de Souza, Ronaldo, 2013, MNRAS, tmp, 2233 
Melbourne, J.; Williams, Benjamin F; Dalcanton, Julianne J.; 
Rosenfield, Rhilip; Girardi, Leo; Marigo, R; Weisz, D.; Dolphin, 
A.; Boyer, Martha L.; Olsen, Knut; Skillman, E.; Seth, Anil C., 
2012, ApJ, 748, 47 

Mouhcine, M. & Lanfon, A., 2002, A&A, 393, 149. 

Moorwood, A. F. M.; Glass, 1. S., 1982, A&A, 115, 84 
Noel, N. E. D.; Greggio, L.; Renzini, A.; Carollo, C. M.; Maras¬ 
ton, C., 2013, ApJ, 772, 58 

Rroctor, R. N., Sansom, A. E., & Reid, 1. N. 2000, MNRAS, 
311,37 

Ramos Almeida, C.; Perez Garcia, A. M.; Acosta-Pulido, J. A., 
2009, ApJ, 694, 1379. 

Ramos Almeida, C.; Perez Garcia,, A. M.; Acosta-Pulido, J. A.; 


Rodriguez Espinosa, J. M.; Barrena, R.; Manchado, A., 2006, 
ApJ, 645, 148 

Rieke, G. H., Lebofsky, M. J., Thompson, R. L, Low, F. J., Toku- 
naga, A. T, 1980, ApJ, 238, 24 

Riffel R., Pastoriza M. G., Rodrguez-Ardila A., Maraston C., 
2007, ApJ, 659, L103 

Riffel, R., Pastoriza, M. G., Rodriguez-Ardila, A. & C. Maraston, 
2008b, MNRAS, 388, 803. 

Riffel, R.; Pastoriza, M. G.; Rodriguez-Ardila, A.; Bonatto, C., 
2009, MNRAS, 400, 273 

Riffel, Rogemar A.; Storchi-Bergmann, Thaisa; Riffel, Rogerio; 
Pastoriza, Miriani G., 2010, ApJ, 713, 469. 

Riffel, R. & Vale, T. B., 2011, Ap&SS,334,35L 

Riffel, Rogfio; Riffel, Rogemar A.; Ferrari, Fabricio; Storchi- 
Bergmann, Thaisa, 201 la,MNRAS, 416, 493. 

Riffel, R.; Ruschel-Dutra, D.; Pastoriza, M. G.; Rodrguez-Ardila, 
A.; Santos, J. F. C., Jr.; Bonatto, C. J.; Ducati, J. R., 2011b, MN¬ 
RAS, 410m 2714. 

Riffel, R.; Bonatto, C.; Cid Fernandes, R.; Pastoriza, M. G.; Bal- 
binot, E., 2011c, MNRAS, 411, 1897 

Rayner, J.T., Cushing, M.C., & Vacca, W.D., 2009, ApJS, 185, 
289 

Salaris, M.; Weiss, A.; Cassara, L. R; Piovan, L.; Chiosi, C., 2014, 
A&A, 565. 

Sanchez-Blazquez, R; Gorgas, J.; Cardiel, N.; Gonzalez, J. J., 
A&A, 457,809. 

Sarzi, M., Rix, H.-W., Shields, J. C., Ho, L. C., Barth, A. 
J.,Rudnick, G., Filippenko, A. V., Sargent, W. L. W. 2005, ApJ, 
628, 169 

Serra, R, Trager, S. C., Oosterloo, T. A.,; Morganti, R. 2008, 
A&A, 483, 57 

SiFchenko, O. K., 2006,ApJ,64L229. 

Storchi-Bergmann, Thaisa; Riffel, Rogemar A.; Riffel, Rogerio; 
Diniz, Marlon R.; Borges Vale, Tiberio; McGregor, Reter J., 
2012, ApJ, 755, 87. 

Tonini, Chiara; Maraston, Claudia; Devriendt, Julien; Thomas, 
Daniel; Silk, Joseph, 2009, MNRAS, 396, 36. 

van der Wei, A.; Franx, M.; Wuyts, S.; van Dokkum, R G.; Huang, 
J.; Rix, H.-W.; Illingworth, G. D., 2006, ApJ, 652, 97 

Zhang, Y; Gu, Q.-S.; Ho, L. C., 2008, A&A, 487, 177. 

Zibetti, Stefano; Gallazzi, Anna; Chariot, Stephane; Rierini, 
Daniele; Rasquali, Anna, 2013, MNRAS, 428, 1479. 


© RAS, MNRAS 000.[T]-?? 


